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DYNAMICS AND DISEQUILIBRIUM CARBON CHEMISTRY IN HOT JUPITER ATMOSPHERES, WITH APPLICATION 

TO HD 209458B 
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ABSTRACT 

Chemical equilibrium considerations suggest that, assuming solar elemental abundances, carbon on HD 
209458b is sequestered primarily as carbon monoxide (CO) and methane (CH4). The relative mole fractions 
of CO(g) and CHi(g) in chemical equilibrium are expected to vary greatly according to variations in local tem- 
perature and pressure. We show, however, that in the p = 1-1000 mbar range, chemical equilibrium does not 
hold. To explore disequilibrium effects, we couple the chemical kinetics of CO and CH4 to a three-dimensional 
numerical model of HD 209458b's atmospheric circulation. These simulations show that vigorous dynamics 
caused by uneven heating of this tidally locked planet homogenize the CO and CH4 concentrations at p < 1 bar, 
even in the presence of lateral temperature variations of ~ 500-1000 K. In the 1-1000 mbar pressure range, 
we find that over 98% of the carbon is in CO. This is true even in cool regions where CH4 is much more 
stable thermodynamically. Our work shows furthermore that planets 300-500 K cooler than HD 209458b can 
also have abundant CO in their upper layers due to disequilibrium effects. We demonstrate several interesting 
observational consequences of these results. 

Subject headings: planets and satellites: general, planets and satellites: individual: HD 209458b, methods: 
numerical, atmospheric effects 



1. INTRODUCTION 



Char bonneau et alJ (2000) and Henrv et al. (2000) discov- 
ered the first known transiting extrasolar giant planet (EGP), 
HD 209458b. Transits are possible if the orbital plane of the 
planet is tilted along our direct line of sight. The detection 
of HD 209458b in transit has been a watershed in the study 
of planets around other stars because physical properties of 
the planet can be inferred directly from the transit data. The 
planet's orbit is nearly circular, with a period of 3.5257 days 
(lLaughlin et alJ2005ah . Furthermore, Laug hlin etal] f2005b) 
have combined transit data with high-resolution spectra of the 
star to infer accurate values for the planet's mass and radius: 
M p = (0.69 ± 0.05)M Jupiter and R p = (1.32 ± 0.05)tf Jupiter . 

Photons from HD209458b and the transiting planet 
TrES-1 have recently been detected dPeming et all [2005b: 
iCharbonneau et ail |2005t). In addition, transmission spec- 
troscopy measurements, in which observers track changes 
in the star's spectral features while the planet is in tran- 
sit, have unveiled key f eatures of the com p osition of HD 
209458b's atmosphere dBrown et all 120011 iBrownl 12001 
ICharbonneauetal] 12001 IVidal-Madiar et alJ 120031 2004). 
The mor e recent detections of transiting planets OGLE- 
TR-56b (iKonacki et alJ 12001 OGLE-TR-113b and 132b 
( iBouchv et al J 12004 TrES - 1 (lAlonso et alJ 12004. OGLE- 
TR-ll lb dPont et alJ l200l . OGLE TR-l Ob (IBouchv et al J 
I2005at IKonacki et al1l2005l). HP 14902 6b (iSato et al]l2005ir 
and HD 189733b (Bouchv et al. 2005b]) reinforce the impor- 
tance of the transit method for the characterization of close-in 
EGP atmospheres. 

HD 209458b resides in close proximity (0.046 AU) to a 
Sun-like star in the constellation Pegasus. The intense stellar 
heating of close-in E GPs drives global-scal e atmo s pheric dy- 
namic s. As shown b y Showman & Guillot ( 2002) . iCho et alJ 
( 120031) . iBurkert et alJ ( I2005I) . and iCooper & Showmanl 
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J20051) . winds of ~ 3000 m s 1 carry material across a hemi- 
sphere in only Td yn ,h ~ 10 5 sec. Vertical transport is rapid as 
well, with typical wind speeds w ~ 20 m s" 1 at p = 1 bar (see 

Theory predicts that HD 209458b is locked by strong stel- 
lar tides to a synchronously rotating state, owing to the short 
time for s pin-orbit sy nchronization of a planet at 0.046 AU 
(see e.g., Barnes & Fortnev 2003). The planet therefore re- 
ceives steHarjadiation on one he misphere only ( iGuillot et alJ 
U996tlShowman & Guillol2002l) . Various effects causing de- 
viations from synchronous rotation have been explored in the 
literature, such as non-zero eccentricity or su bstantial obliq- 
uity pumped by the presence o f another planet ( Laughli n et alJ 
2005b; Winn & Holman 2005). Howe ver, the apparent lack o f 
a third massive body in the system (Laughlin et al. 2005a), 
as well as the t iming of the transit and secondary eclipse 
dPeming et all2 005b). makes a large eccentricity or obliquity 
unlikely for this system. 

We study the atmosphere under the assumption of syn- 
chronous rotation. We note, however, that deviations from 
a completely synchronous state are conceivable. For exam- 
ple, effects of atmospheric tides a re important on Venus, as 
shown by Correia & Laskai (2001). In an order-of-magnitude 
calculation, Showman & Guillot (2002) suggested that atmo- 
spheric dynamics could induce deviation from synchronous 
rotation by up to a factor of 2. A factor of 2 is the extreme 
case; 10-20% deviation might be more typical. This estimate 
depends on the planet's tidal dissipation factor (0, as well 
as properties of the circulation itself and the processes trans- 
ferring angular momentum and rotational kinetic energy from 
the atmosphere to the interior, which are highly uncertain. 

The GOV primary irradiates the planet up to an extremely 
high effective temperature (T e ff ~ 1340 K; see llro et alJ2005t) 
compared to the planets of our Solar System. This leads to 
very different chemistry than we are familiar with in plan- 
etary atmospheres. The atmosphere is primarily composed 
of hydrogen and helium. Carbon, nitrogen, and oxygen 
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are all minor constituents. Despite their low abundances, 
these elements are important for determining the spectral 
characteristics of the planet beca use they form molecule s 
having strong spectral signatures (Seager & Sasselov 1998). 
Radiative-equilibrium calculations bv llro et all {2005) show 
that in the 1 mbar to 10 bar pressure regime, which is the 
region of primary interest for observations, the temperature 
lies in the 1000-2000 K range. At these temperatures and 
pressures, chemical equilibrium considerations suggest that 
carbon on HD 209458b is sequeste red primarily in c arbon 
monoxide (CO) and met hane (CH4) dBurrows & Share! 19991: 
Lodders & Feglev 2002). These models show that CO(g) is 
the most stable carbon-bearing species at high temperatures 
and low pressures. At low temperatures and high pressures, 
CH4(g) is more stable thermodynamically. The oxygen not 
bound up in rock is present as either CO(g) or water vapor, 
H 2 0(g). 

In planetary atmospheres, however, chemical equilibrium 
may not hold. There are several reasons for this. First, photo- 
chemical dissociation can elevate the concentrati ons of pho- 
toche m ical products above eq uilibrium levels ( Mo ses et alJ 
2000). iPrinn & Barshavl i ll 9771) note that photochemical ox- 
idation of CH4 is a well-known source of CO in the Earth's 
atmosphere. Second, impacts by meteorites, comets, etc., de- 
posit metals and oth er compounds into the atmospher es of gi- 
ant planets (see, e.g., Pra ther et alJl 978: Moses 1992). Third, 
chemical reactions may not always occur quickly enough rel- 
ative to dynamical mixin g for chemical equ ilibrium to be at- 
ta ined JPrinn & Barshavll977HSmithll998l) . 

Liang et al. ( 2004) show that the first effect — 
photochemical production — is unlikely to significantly 
alter the concentration of CO in HD 209458b's stratosphere 
(p < 1 mbar). It is unknown how much impactors have 
modified the composition of HD 209458b's atmosphere 
durin g its ev o lution . We note, however, that calculations by 
iBezard et alJ (|2002) show that CO deposition by impactors 
cannot explain the observed abundance of CO (~ 1 ppb, as 
opposed to essentially zero in equilibrium) in Jupiter's lower 
troposphere. On Jupiter, rather, disequilibrium results from 
convective upwelling from deep regions (p ~ 300 bar), where 
the equilibrium abundance of CO is higher. As the reaction 
rates plummet with decreasing temperature and pressure in 
ascending fluid parcels, the CO concentration quenches to its 
observed value. 

Photospheric pre ssures lie in the 10-500 mbar range for 
all IR wavelen gths (Fortnev et alJl2005t lBurrows et alJl2005t 
iBarman et alJ 1200*51 Seageretal. 2005). If large tempera- 
ture variations of 500-1000 K in th is reg ion are a reality 
on HD 209458b dShowman & Guillotll2002t ICho et alJl2003l 
iBurkert et alJl2005t TCooper & Showman 2005), it is conceiv- 
able that chemical quenching occurs in the horizontal rather 
than the vertical direction. In this hypothetical scenario, fluid 
parcels attain chemical equilibrium on the planet's dayside, 
where the temperature is high and the reaction kinetics are 
fast. Vigorous dynamical motions then carry these equili- 
brated parcels to the nightside, where they rapidly cool to 
temperatures where the reactions are too slow to equilibrate 
the system. Depending on the meteorology and variations of 
the chemical rates with temperature and pressure, complex 
lateral CO distributions could potentially result. For example, 
it would then be possible for CO to be abundant at the evening 
limb but depleted at the morning limb. 

To test the hypothesis that disequilibrium processes modify 
the abundances of CO and CH4 in the observable atmosphere 



of HD209458b, we couple the reaction kinetics to a three- 
dimensional nonlinear model of the atmospheric dynamics. 
We show that at the temperatures and pressures relevant for 
observations of HD 209458b's atmosphere, the timescale for 
CO/CH4 interconversion is everywhere long: T c h em ^S> 10 5 s. 
Disequilibrium results in regions where the timescale for ad- 
vection of CO by dynamics, Td yn , is shorter than r c h e m- Due to 
uniformly large values of T C h em for p < 1 bar, the simulations 
definitively show that quenching occurs in the vertical, not the 
horizontal. As on Jupiter, we find that the concentration of CO 
for p < 1 bar depends on its value deeper in the atmosphere. 
Vertical quenching and efficient horizontal mixing lead to a 
high-abundance, homogeneous distribution of CO at the pho- 
tosphere. As a result, in cool regions (e.g., the nightside), the 
CO concentration exceeds equilibrium values by many orders 
of magnitude. 

Curiously, CO was not detected on HD 209458b in a re- 
cent g round-based se arch using transmission spectroscopy at 
2/xm llDeming et alJl200"5ah . Several possible reasons for the 
CO null detection have been noted by the authors and others. 
First, CO may not be abundant in this planet's atmosphere, 
though this is only possible if the planet is significantly cooler 
(by ~ 5 00 K) than predicted by current radiative-e quilibrium 
mode ls ilBurrows et al J 120041: IChabrier et~al]l2004t llro et alJ 
120051) . Second, CO may present in th e visible atmosp here in 
some regions but depleted at the limb tiro et al. 2005). Third, 
the spectral signature of CO at the planet's limb may be ob- 
scured eithe r due to photochemical haze or a high-altit ude sil- 
icate cloud (Fortnev et ayKQS^LISyfflBi^SllESQa)- This 
explanation is favored by Deming et al. (2005a). It should 
be noted, however, that the CO null detection reported by 
Deming et al. (2005a) relies on the radiative-equilibrium cal- 
culations of Sudars kv et alJ ll2003l) . If unexpected opacity 
sources (though not necessarily clouds) play a strong role 
in the formation of spectral lines near 2 /im, an upper limit 
to the CO abundance cannot be made with high confidence 
based on these observations. We note CH4 has not been seen 
on this planet either, desp ite attempts to detect its signature 
( iRichardson et all2003b1al) . 

We review in 5j2]the theoretical reasons to expect CO in high 
abundance in HD 209458b's atmosphere based on chemical 
equilibrium arguments. We discuss in 5j5]the main chemical 
processes involved and derive the timescale for CO/CH4 in- 
terconversion. In &\ we describe our coupl ing of chemistry 
to the atmospheric dynamics model of ICooper & Sho wman 
( 120051) . We present the results of these simulations in |5] We 
explain and interpret these results in |6j comparing o ur sim- 
ulation s to a simpler view of quenching presented by Smith 
( 1998). We enumerate the major uncertainties in these calcu- 
lations in §Q In HI we discuss implications of this work for 
near-future observations of HD 209458b and other close-in 
EGPs. 

2. CHEMICAL EQUILIBRIUM CONSIDERATIONS 

Prior discussions of chemistr y on EGPs are based 
on chemical equilibrium models ( iBurrows & Shard 119991: 
Lodders & Feglev 2002). Solar abundances of the elements 
(Anders & Grevesse 19891 lLodders & Feglevll 19981 ) are gen- 
erally assumed in these schemes, which perform global mini- 
mization of the Gibbs free energy of the system, often includ- 
ing thermochemical data for hundreds of compounds. This 
minimization algorithm allows the concentrations of all chem- 
icals to be solved for simultaneously. In such thermochemi- 
cal calculations, once the elemental abundances in the atmo- 
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sphere have been specified, the abundances of all chemical 
species in the system depend solely on the local temperature 
and pressure (allowing f or simple ideas li ke rain out). 

Simulations by Cooper & Showman (2005) suggest that 
the temperature at the photosphere of HD 209458b varies 
across the globe from ~ 1000-1500 K. Chemical equilibrium 
models show that at these temperatures, the atmosphere of 
this planet sho uld contain app reciable quantities of CO and 
CH 4 . lLodders & Fedevl 12002) show that HD 209458b's at- 
mosphere straddles the line of equal CO/CH4 abundance, in 
the absence of disequilibrium effects. The relative concentra- 
tions of CO and CH4 are strong functions of local temperature 
and pressure. At 1 bar, where T < 1 150 K, CH4 is the domi- 
nant carbon-bearing species (in chemical equilibrium), while 
CO is more abundant at this pressure for higher temperatures 
(T>1150K). 

In the temperature and pressure ranges of interest for this in- 
vestigation, calculation of the equilibrium abundances of CO 
and CH4 is straightforward and reduces to several simple alge- 
braic expressions. These expressions are used in our disequi- 
librium model, in which the mole fract ion o f CO is relaxed to 
its value in chemical equilibrium (see ^4.2> . The expressions, 
as we will show, include the metallicity dependence directly. 
This allows us to explore the sensitivity of our disequilibrium 
calculations to the metallicity of HD 209458b's atmosphere, 
which is conceivably enhanced from that of the host star. Such 
enhancement is seen on Jupiter, where C, N, and probably O 
are several times more abu ndant relative to hy drogen than at 
the photosphere of the Sun ( Atrev aet alll999l) . 

We define X[ to be the mole fraction of the i-th species in 
the system; e.g., Xn 2 is the mole fraction of molecular hydro- 
gen. It is given by the ratio of the partial pressure of the i-th 
species to the total atmospheric pressure (i.e., Xi = pi/ p). At 
the temperatures and pressures in this region of HD 209458b's 
atmosphere, hydrogen is primarily molecular H2 [J. Fortney, 
private communication]. We take X^ = 0.83 for the H2 mole 
fraction throughout this paper (i.e., we assume that H2 disso- 
ciation is everywhere negligible). 

Although oxygen is mostly in CO and H2O, some oxygen 
is al so sequestered in rock, most of w hich forms deeper than 1 
bar (Lodders & Feglev 2002; Fortnev et al. 2005). We assume 
16% of the atmospheric oxygen (the maximum possible in a 
solar abundance mixture is in silicates. Depending on where 
in the atmosphere MgSi03 condenses, this may underestimate 
(by ~5%) the oxygen available to affect CO/CH4 chemistry 
in the region of interest. 

Defining ci and C2 as the mixing ratios of gaseous carbon 
and oxygen in the atmosphere, the mole fractions Xco, Xqh 4 , 
and XiinO are related through mass balance: 



X C h 4 +Xco = ci=5.91 x 10" 4 
Xco+X H2 o = c 2 = 1.04x 10" 3 



(1) 

(2) 



for solar metallicity. We see from equations[2and[2]that wa- 
ter vapor will always be present in the system, even if 100% 
of the carbon is in CO, since oxygen is more abundant than 
carbon. 

The third expression relating the quantities Xqo, Xch 4 , and 
Xh 2 o derives from the equilibrium expression of the net ther- 
mochemical reaction 

CO + 3H 2 <-> CH4 + H2O. (3) 

The equilibrium expression for reaction [3] which we denote 



= exp 



AfGco - A t GcH 4 _ AfGn 2 o 
RT 



jeq 



XcoX^p 2 

Xch 4 Xh 2 

(4) 

In equation |3 p is atmospheric pressure [bar], T is the 
temperature [K], and R = 8.314 J mol" 1 K" 1 is the universal 
gas constant. The symbols AfGco, AfGcH 4 , and AfGH,o 
represent the Gibbs free energies of formation of CO, CH4, 
and H2O, respectively. The Gibbs free energies of formation 
depend on temperature sensitively but only weakly on pres- 
sure. We simply use the values at standard pressure (1 bar), 
which are tabulated in the NIST-JANAF Thermochemical Ta- 
bles dChasell998l) . 

With the three equations \l\ [2] and [4] we can solve for the 
three equilibrium mole fractions Xco, Xu 2 o, and Xch 4 as a 
function of atmospheric pressure (p) and temperature (T). 
Defining / = f^ eq /(p 2 X^), we sorve f° r ^co- We take the 
negative root of the quadratic equation, which represents the 
physical solution: 



-b—\/b 2 — 4ac 
2a 



(5) 



with a = f, b = —f(c\ +cz)— 1, and c = fc\C2- The other 
concentrations are then easily solved for using the mass bal- 
ance expressions (equations \l\ and 0: Xch 4 = c\ - Xco and 
^h 2 o = C2~Xco- The solution is valid for temperatures be- 
tween 500-2500 K over a broad range of pressures in the at- 
mosphere. It is also valid for atmospheres with enhanced or 
depleted C and O abundances by simply tuning c\ and C2 (so 
long as c\ and C2 are kept small: ci, C2 <C 0.1). 

We show in Table the ratio of CO to total carbon in the 
atmosphere (Xco / ci), assuming chemical equilibrium, for 
a variety of temperatures and pressures. The table shows 
that the CO abundance — assuming equilibrium conditions — 
is a strong function of pressure and especially temperature. 
If large temperature variations in the atmosphere are indeed 
a reality on HD 209458 b, as pre dicted both in radiative- 
equilibrium models (e.g., |IroetalJ 2005) and si mulat ions of 
the atmospheric dynamics ( Cooper & Showman 2005), large 
gradients in the equilibrium concentrations of CO and CH4 
result. We now consider disequilibrium effects, which lead to 
vastly different conclusions about the relative concentrations 
of CO and CH4 in the upper layers of this planet's atmosphere. 

3. CHEMICAL PATHWAYS FOR CO REDUCTION 

CO concentrations above equilibriu m valu e s have been 
disc overed both on Jupiter ( iBeeJ 1 19751 lLarson et all 
19781) and on the brow n dwar f Gl 229B dNolletal] 
1997t lOppenheimer etall 119981 iSaumon et all 12 000). 
Feglev & Lodders (1996) and iGriffith & Yelld J 19991) show 
that CO disequilibrium in Gl 229B results from the inef- 
ficiency of chemical reduction of CO in regions of low 
temperature and pressure relative to the rate of vertical 
mixing. The CO seen in the spectrum of Gl 229B originates 
from deeper in the atmosphere, where CO is abundant in 
chemical eq uilibrium. 

Similarly. Iprinn & BarshavH 19771) suggest that the CO de- 
tected in the troposphere of Jupiter was transporte d up by con- 
vectio n from deeper, hotter layers of the planet. Beza rd et alJ 
(2002) explore several explanations for the unexpected 1 ppb 
CO concentration observed on this planet. They conclude that 
convective mixing from deeper layers remains the most likely 
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source of Jupiter's tropospheri c CO, although the reactions 
suggested by Prinn & Barshav ill 9771) are now thought to be 
too slow to explain the observed CO abundance. 

We hypothesize that dynamics may drive the abundance of 
CO on HD 209458b away from chemical equilibrium values. 
To ascertain the extent to which dynamics perturbs chem- 
ical equilibrium in this atmosphere, we first calculate the 
timescale for the processes converting CO to CH4 and back. 
This depends on the chemical reactions involved. Two inde- 
pendent reaction pathways for the chemical reduction of CO 
have been considered in the literature. 

The following reaction was suggested bv lPrinn & Barshav 

H 2 + H 2 CO — * OH + CH 3 (6) 

This reaction is the rate-limiting step in a three-part sequence 
resulting in the chemical reduction of CO to CH4. These au- 
thors calculate a rate constant for reaction[6]of 



te]=2.3 x 10" 10 exp 



36, 200 K 



cm 3 s 1 . 



(7) 



tPrinn & Barshav! i 19771) note that the other steps leading to 
the reduction of CO to CH4 are fast relative to reaction [6] 
They use equation[7]to estimate the chemical lifetime of CO 
in disequilibrium as a function of temperature and pressure. 

Based on more recent measurements, however, 
iGriffith & Yelld i 19991) discuss that the reaction rate for 
reaction \Q is sig nificantly slower than the rate assumed by 
iPrinn & B arshaj] (119771) : 



:3.4x 10" s - 



-1.12 



exp 



43,192K 



cm 3 s 



(8) 

convective 
derive the 



Furthermore, as shown by Smith (1998), the 
timescale used by Prinn & Barshav ( 1977) to 
quench level is incorrect because it does not account for the 
gradient in both the lifetime and equilibrium abundance of CO 
over an atmospheric scale height. Using the updated reaction 
rate, equationlHl and the correct timescale for convective mix- 
ing, IBezard et a find that the reaction sequence for 
CO destruction suggested bv lPrinn & Barsha^ ( 119771) cannot 
explain the disequilibrium abundance of CO in Jupiter's tro- 
posphere. A faster reaction pathway for the reduction of CO 
to CH4 must be soug ht. 

Yun g et aT] (1198 8) investigate an alternative three-step pro- 
cess involving the methoxy radical CH3O. They argue this 
pathway is more kinetically favorable because the strong car- 
bonyl (C = O) bond in H 2 CO is broken in a separate reaction 
step from the form ation of the four C-H bonds in methane. 
lYung etaD (fT988l) propose the following three-stage reaction 
set: 

H + H 2 CO + M -> CH3O + M (9) 



CH3O + H 2 -> CH3OH + H 
CH3OH + H -> CH 3 + H 2 



(10) 
(11) 



In this sequence, the first reaction, in which the C = O bond 
is broken, is the rate-limiting step. The net reaction resulting 
from summing reactions l9l 1 101 and ll 1 l is 



H 2 CO + H 2 + H -> CH 3 + H 2 0. 



(12) 



We derive the number density of H 2 CO by assuming equi- 
librium with H 2 and CO: H 2 + CO <-> H 2 CO, which is fast 
compared to reaction |9] We note also that the production of 



methane (CH4) from CH3 occurs with nearly 100% efficiency 
JYung et alJ1988l) . 

The rate of reaction|9]has never been measured i n the labo- 
ratory , but the rate of the reverse reaction is known ( Pa ge et afl 



1 R = L4X1 °" 6 VK 



-1.2 



7,800K\ 
exp [ — — ) cm' s . (13) 



EauationflJlrepresents the low-pressure limit of the rate of the 
reverse of reaction|5]in the temperature range relevant for this 
plan et. For the high-pre ssure limit, we adopt equation (11) 
from lBezard et alJ (12002 ) : 



= 1.5 x 10 1 



■ exp 



12,880K 



(14) 



The net rate constant for the reverse reaction |9jR can be ex 
pressed as 



(15) 



where n [cm 3 s i here is the n umber density of molecules in 
the atmosphere (Bezard et al. 2002). It is related to pressure 
through the ideal gas equation of state. 

Using the proper eddy mixing timescale for chemical 
que nching derive d bv ISmithl |l998), Griffit h & Yeliel (119991) 
and Bezar det alJ (120021) demonstrat e that the kinetic rates 
derived from the Yung et al. (1988) reactions lit observa- 
tions of CO in disequilibrium on the brown dwarf Gl 229b 
and on Jupiter. As not ed above, the reaction proposed by 
IPrinn & Barshay ill 9771) is too slow to account for the ob- 
served cojic^ntratiojis ofCO in t hese objects. We th erefore 
follow lGriffith & Yelld i ll 9991) and lBezard etail ll2002l) in this 
treatment and assume that the catalyzed Yung et al. ( 1988) re- 
actionsl9lll Hare the most efficient mechanism for the conver- 
sion of CO to CH 4 . 

For the purpose of determining the CO quench level on 
HD 209458b, we use the measured reaction rates (eq. 1131 - 
1 1 51 along with thermochemical data to determine r c h e m- This 
quantity, which is a strong function of both pressure and tem- 
perature, represent s the timescale over w hich CO can remain 
in disequilibrium (Prinn & Barshavl ll977l) . Conservation of 
mass (eq. Q dictates that the rate of CO destruction must 
equal the rate of CH4 production. So, r C hem also represents the 
timescale for conversion of CH4 back to CO. This would not 
be so if other molecules competed effectively for the available 
carbon atoms. At these pressures and temperatures, though, 
CO and CH4 are the predominant carbon-bearing species (see 
©. Hereafter, we simply refer to r C hem as the CO/CH4 inter- 
conversion timescale. 

The CO/CH4 i nterconversion timescale is given by 
IBezard et all I2C 



7"chem 



Xco 



-dXco/dt n 2 Am- X H X HlC o ' 



(16) 



where ^jis the reaction rate of|5] An extra factor of n arises in 
the denominator from writing the expression in terms of CO, 
H, and H 2 CO mole fractions rather than concentrations. 

As has been noted, however, only the rate of the reverse of 
reaction[9]has been measured. We therefore use the equilib- 
rium constant of the chemical reaction 



CH3O <-> CO + 3/2H 2 



(17) 
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to express r c h em in terms of equation^] The equilibrium con- 
stant of reactionfTTlis 



eq 



= exp 



A f G(CO)-A f G(CH 3 0) 
RT 



(18) 



The Gibbs free energy of formation of CH3O cannot be found 
in the JANAF tables. W e use values for A f G(CH 3 0) from 
iTsang & Hampsonllll986l) . 

Using these thermoc hemic al data along with equation IT~5l 
we can express r c h e m in a form that can be evaluated numeri- 
cally: 



7"chem(P,T) : 



Xco 



■k h T 



«^R- X CH 3 /|9] r Xh{V/2' 



(19) 



where k\, is Boltzmann's constant. The first equality in equa- 
tion^] is analogous to equation [TBI but for reaction |9] in the 
reverse direction. All intermediate reactions in the system are 
assumed to be in equilibrium. The second equality results 
from substituting equationlT8"1in for the ratio ^co/^ch 3 o- 

We show values for T c h em in Table |2] at the same pressures 
and temperatures chosen in Table Table |2] demonstrates 
that for p < 1 bar, r chem > 10 5 s for all T < 2000 K, which 
is an upper limit to the temperature at HD 209458b's photo- 
sphere. This result immediately rules out horizontal quench- 
ing as an important disequilibrium process at the photosphere. 
The timescale for horizontal mixing Td yn .h ~ 10 s, so that 
Tchem 3> Td yn .h on both the day and night sides of the planet. 

Deeper down (p > 10 bar), CO/CH4 interconversion is fast 
(Tchem <C 10 5 s). At this pressure, the chemical interconversion 
timescale is less than both the horizontal and vertical mixing 
times. We would therefore expect chemical equilibrium of 
CO/CH4 to hold in the atmosphere at p > 10 bar. The top of 
this region (p ~ 10 bar) is the approximate "quench level." 
Note that the precise altitude of quenching also depends on 
atmospheric temperat ure and wind speeds, w hich are them- 
selves highly variable (Cooper & Showman 2005). As a con- 
sequence, the quench level is not an exact isobar. Referencing 
Tabled we see that abundant CO at the photosphere of this 
planet is likely because vertical tran sport processes are prob- 
ably efficient (IShowman & Guillotl 120021) . and the chemical 
equilibrium CO concentration at the quench level is likely to 
be quite high. 

We now discuss a coupled atmospheric dynamics and dis- 
equilibrium CO/CH4 chemistry model to predict CO concen- 
trations for a variety of atmospheric conditions. Our simula- 
tion results (see ^5} verify the general validity of these con- 
siderations. 

4. COUPLED CHEMISTRY/DYNAMICS MODEL 

4. 1 . The Atmospheric Dynamics Model 

We adapt the model of ICooper & Showmanl (120051) to in- 
clude simple carbon chemistry based on the discussions in 
33 and 50 Our model emplo ys the ARIES/GEOS D ynam- 
ical Core, version 2 (AGDC2:" lSuarez & Takacsll 19951) . The 
AGDC2 solves the primitive equations of dynamical mete- 
orology, which are the foundation of n umerous climate and 
nume rical weather prediction models ( Holto*ni ri992t iKalnavl 
2003). 

The primitive equations are a simplification of the full 
Navier-Stokes equations of fluid mechanics, which assume 
hydrostatic balance of each vertical column of atmosphere. 
This is equivalent to assuming that the vertical acceleration 



is negligible compared to the buoyancy. The approximation 
holds in stably stratified regions for "shallow" flow, in which 
the vertical de pth is much smaller than the horizontal extent 
( Kalnav 2003, Ch. 2). Heavy irradiation from the star extends 
the radiative zone of this pla net down to p ~ 1 kbar. This 
is ~1/10 of the radi us (iBurrows et al] 120031 IChabrier et alJ 
2004; Iro et al. 2005), leading to a horizontal to vertical as- 
pect ratio of about 10: 1 . 

The AGDC2 is a finite-difference model, which discretizes 
the d ynamical variables onto a staggered latitude-longitude C- 
grid ( Arakawa & Lamb 1977). The AGDC2 consists of a set 
of subroutines that compute the time tendencies of the prog- 
nostic variables, inclu ding the zon al and meridional winds, 
potential temperature (Holton 1992, p. 52), surface pressure, 
and an arbitrary number of passive tracers. At each time step, 
the AGDC2 simply updates the time tendencies to include the 
effects of the dynamics. All physical tendencies (e.g., chem- 
istry, radiation, etc.), as well as temporal and spatial filtering 
needed to control computational instabilities, are performed 
outside of the AGDC2. 

It is not possible in general to find numerically sta- 
ble solutions to the primitive equations without viscosity 
dDowling et al.l H998I) . The standard approach, which we 
adopt here, is to add diffusive terms to the momentum, en- 
ergy conservation, and tracers equations. We have experi- 
mented with the diffu s ion pa rameters, espousing the philos- 
ophy of Polvan iet alJ d2004t) to use as weak hyperdiffusion 
as possible consistent with numerically stable solutions. For 
the simulations discussed here, we use 4 th -order hyperdiffu- 
sion, with coefficients set to damp grid-scale instabilities over 
a timescale of ~ 30 minutes. 

For all simulations, we take the "top" layer of the model to 
be at 1 mbar. Although a proper model of the upper atmo- 
sphere would require the model top be placed at a lower pres- 
sure, we generally ignore the upper atmosphere in this work. 
Modeling the high atmosphere would require a proper treat- 
ment of stratospheric photochemistry and the possible effects 
of non-local thermodynamic equili brium in the radiation field 
at su b-millibar pressures (Liang et al. 2004; Barm anet alJ 
2002). The atmosphere spans ~ 15 pressure scale heights be- 
tween our input top pressure (1 mb) and the planet's radiative- 
convective boundary (1 kbar). By comparison, the pressure in 
the Earth's atmosphere spans only 2-3 pressure scale heights 
from the surface to the tropopause. We use 40 layers evenly 
spaced in log pressure. 

We assume in this treatment that synchronous rotation holds 
in the convective interior. The consequences of this assump- 
tion are twofold: (1) the rotational period is fixed to the 
3.5 day orbital period (leading to a moderately strong Cori- 
olis force), and (2) heating by the star occurs on one hemi- 
sphere of the planet continually. This assumption is valid 
for many hot Jupiters (such as HD 209458b), which expe- 
rienc e rapid tidal circularization of their orbits (Guillot et al. 
1996). Atmospheric tides may cause the orbit to be sl ightly 
asynchronous, as noted by Showman & Guillot (2002). But 
large deviations from synchronous rotation are unlikely on 
HD 209458b without t he prese nce of another (thu s far unseen) 
planet (Laughlin et al. 2005a). We discuss in 37.21 that the 
asynchronous component (a factor of 2 in the extreme case, 
though more likely in the 10-20% range) will not greatly af- 
fect the outcome of our predictions for CO/CH4 chemistry 
(see {J7J. But it would be interesting in future studies to ex- 
plore the dynamics of hot Jupiters under the alternate assump- 
tion of variable heating due to non-zero eccentricity, large 
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obliquity, or strong the rmal tides. 

We follow ICooper & Showmanl d2005t) for the scalar input 
parameters, ignoring small variations in the gravitational ac- 
celeration, heat capacity, and mean molecular weight in the 
atmosphere. We set these values to constants: g = 9.42 m s~ 2 , 
[1 = 1 .8 1 x 10" 3 kg moF 1 , and c p = 1 .43 x 10 4 J kg" 1 K" 1 . The 
values of these constants have been chosen to be intermedi- 
ate between the values they would have at 1 mbar and 3 kbar 
because H2 dissociation is expected to occur in the deepest 
layers of the model. For the chemistry, we assume hydrogen 
is molecular everywhere, which holds in the upper layers of 
interest here (see S0. 

4.2. Coupling CO /CH4 Chemistry 

The AGDC2 includes horizontal advection for an arbitrary 
number of passive tracers. We denote these as q^, with 
k = 1,2, ...,km. As few as possible is desirable, as the tracer 
equations add appreciably to the total computations required 
per time step. 

Tracers are passive in the sense that they are advected 
by the flow but do not themselves influence the dynamics. 
Unless interdependencies are added, the tracer solutions in 
the AGDC2 do not affect the other prognostic variables (or 
each other). The equations operating on each qu are of the 
same form as those used for potential temperature advection 
ISuarez & Takacslll995l) . The advection scheme is fourth- 
order in this version of the dynamical core. 

We add CO as a passive tracer to the AGDC2. The other 
carbon and oxygen bearing species, including CH4, H2O, 
CH3O, HCO, etc., are either very short-lived or indirectly 
obtainable through mass balance (eq. ^ an d|3. To confirm 
thi s, we calculated li fetimes for the intermediate species in 
thelYun g et aTl f 1988) sequence. These were all much shorter 
than our T c h em , which describes the timescale for interconver- 
sion of CO and CH4. For the purposes of this investigation, it 
is sufficient therefore to cons ider just a sing le tracer. 

We follow equation (2) of Smith ( 1998) to express the ef- 
fect of the dynamics on CO/CH4 chemistry in terms of the 
quantities we know, Xco.eq and r c h em - The idea is to relax Xco 
to its value in chemical equilibrium, Xco^a,'- 

DXco _ Xco~ Xco. eq »20) 

Dt 7"chem 

We write D/Dt (not d/dt or d/dt) to emphasize we mean 
the total (or material) derivative, not the local derivative. As 
usual, the advection terms must be included w hen th e equa- 
tions are written in the Eulerian frame llHoltonll 19921 p. 28- 
31). 

In equation |20| we see that if T c h em is short relative to dy- 
namical timescales, chemical equilibrium of CO will be at- 
tained in the atmosphere. In regions where T c h em is longer 
than Tdyn, CO will not be in equilibrium. In the former case, 
the CO abundance will depend solely on the local tempera- 
ture and pressure (which themselves vary greatly throughout 
the region of integration). In the latter situation, whether the 
actual value of Xco is greater than or less than the equilibrium 
value depends on global atmospheric conditions. 

In the troposphere of Jupiter, for example, Xco is extremely 
high relative to Xco.eq because convection mixes CO upwards 
from the quench level (at several hundred bar s), where the 
equili brium abundance of CO is appreciable (Beza rdet al.l 
2002). The situation of HD 209458b is different; the atmo- 
sphere in the 10 mbar to 10 bar region is stable to turbulent 
convection due to the huge flux of stellar radiation from above 



( iGuillot etalJ fl996). As discussed bv IShowman & GuilioH 
(2002), however, vertical velocities of ^20 m s" 1 are plau- 
sible on this planet, implying vigorous vertical transport. The 
CO concentration for p < 10 bar depends critically on the ver- 
tical flow, which we will discuss in 5j6] 

4.3. Newtonian Heating Scheme 

We defer to future work the ambitious task of coupling at- 
mospheric dynamics to a true radiative transfer code. Rather, 
we approximate the effects of radiation using a simple Newto- 
nian heating prescription. In this scheme, the thermodynamic 
heating rate at each grid point is proportional to the difference 
between the prescribed 3D radiative-equilibrium temperature 
and the t rue local temperature . 

As in Coope r & Showmanl J2005I), we start with the ID 
radiative-equilibrium T(p) profile of llro et all (|2005). Their 
radiative-equilibrium profile, which we call Ti m (p), assumes 
heat redistribution over the entire globe. The calculation in- 
cludes the major sources of gaseous opacity: Rayleigh scat- 
tering, collision induced absorption, bound-free and free- 
free absorption, molecular rovibrational bands from H2O, 
CO, CH 4, and T i O, an d resonance lines from the alkali 
metals, llro et alJ (2005) calculated a line-by-line solution 
to the monochroma ti c radi ative transfer equation using the 
Goukenleuaue et al. (2000) model. This code makes the 
two- s tream approximation (see e.g., [Chamberlai n & Hunted 

QUI, 

and the atmosphere is t aken to be plane- parallel. We 
also use values published by llro et al.l 112005) for T ra d(p), 
which is the radiative relaxation time constant. This is an in- 
creasing function of pressure. The value of r ra d is ~ 3 hours 
at 1 mbar (which is the topmost layer of our model) and in- 
creases to ~ 1 year near 10 bar. 

We use the algorithm described in Coo per & Showmanl 
(2005) to calculate T n [ g h t (p) and T ss (p) from T bo (p). These 
are the radiative-equilibrium profiles on the nightside, which 
receives no direct heating, and at the substellar point, where 
the star is directly overhead. We assume in these calcu- 
lations that nightside radiative-equilibrium temperatures de- 
pend only on pressure, not longitude and latitude. In this 
procedure, we prescribe a temperature difference in radiative 
equilibrium between the substellar point and the nightside, 
AT eq (p) = T ss {p)-T„i g j it {p). This is an adjustable parameter, 
which governs the strength of the radiative forcing. Its has 
been chosen to be logarithmic in pressure, so the profiles con- 
verge deep in the atmosphere. 

The substellar and nightside radiative-equilibrium profiles 
are then computed by balancing the total flux radiating to 
space at the top of the atmosphere. The nightside and inte- 
grated dayside contributions must sum to the known globally 
averaged value, AuR^aT^: 

4irR 2 p aT* = 2nR 2 p aT^ t +R 2 p f f ar d 4 ay (A, 4>) cos <^dA. 

J J day 

(21) 

The coordinates (A, <fi) are longitude and latitude, respectively. 

The flux balance in equation |^ yields a transcendental 
equation for r n ; g ht at the top of the atmosphere in terms of 
the free parameter, Ar eq : 

47^ = 37; 4 ght + (r night + Ar eq ) 4 . (22) 

In equation|22] 2 factors of r n 4 oht come from the nightside; the 
integral over the dayside contributes 1 more factor of T^ hv as 
well as the term involving Ar eq . Equation 1221 can be solved 
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FIG. 1 . — Radiative-eq uilibr ium T(p) profiles used in the Newtonian heat- 
ing scheme described in SI4.3I The profi les shown are th e nightside radiative- 
equilibrium profile (long-dashed linel. Ilro et alJ l2005Vs globally averaged 
radiative-equilibrium profile (solid line), and the radiative-equilibrium pro- 
file at the substellar point (short-dashed line). At the topmost layer of the 
model atmosphere (1 mbar), we have set the temperature difference between 
the substellar point and the nightside to 1000 K. The nightside and substellar 
profiles are balanced so the integrated outgoing flux at the top of the atmo- 
sphere equals 47rRpCr7j* . The profiles converge logarithmically with increas- 
ing pressure (depth). 



numerically; we use Newton-Raphson fsee lPress et all 19921) . 
r ni „ ht = 500 K is the solution for Ar eq = 1000 K and 7^ = 
1000 K. 

The above radiation balance is only computed once — at 
the top model layer — to ensure the outgoing fluxes are self- 
consistent. We emphasize this is not a rigorous t r eatme nt 
of radiation. The goal is simply to use Ilro et a D l|2005)'s 
radiative-equilibrium profile, which is one-dimensional, to 
prescribe the Newtonian heating rate over the model grid in 
a way that does not significantly alter the radiation budget of 
the atmosphere. 

We show in Figure^the radiative-equilibrium profiles gen- 
erated by this method from pressures of 1 mbar to 10 bars. 
We plot Tj m (p) as well for comparison. The profiles plotted 
are for our nominal HD 209458b model, in which we have 
set Ar eq = r ss - r n j g ht = lOOO K at the top of the atmosphere 
(see 5j5]for a more complete description). We only show the 
profiles down to the 10 bar level because the radiative time 
constants deeper than that are very long (t^h ^> 10 7 s; see 
Ilro et alJ (12005)). As in Coop er & Showmad 120051) . we have 
simply turned off heating for all pressures greater than 10 bar. 

4.4. Limb Effects 

Two limb effects may be important at low pressures. 
We have not attempted here to model them in detail. 
But to account for them , we use a modified form of 
ICooper & Showman! J2005), equation(2). 

First, near the limb, the angle of incidence of the imping- 
ing radiation is large. The optical path length is longer in the 
slant geometry than at normal incidence. Hence, the stellar 
radiation will not penetrate as deeply near the limb but will 
be absorbed higher up in the atmosphere. This causes addi- 
tional heating at pressures significantly less than 1 bar. Sec- 
ond, the star is not a point source; it has finite angular size 
in the planet's sky. This creates a small region of continuous 
heating past the day-night terminator, which we define here 
as 90° from the substellar point. Here, the disk of the star is 




20 40 



60 



80 100 120 140 160 180 
a 



FIG. 2. — The angular dependence of the 3D radiative-equilibrium profile, 
which is symmetric about the star-planet line. We define a to be the angle 
between the local normal vector at a particular longitude and latitude (A, <f>) 
and a ray from the planet's center passing through the substellar point (A = 0° , 

= 0°). This is a s lightly warmer profile near the limb than that used by 

1 Cooper & Showman 2005). 

partly visible. 

These considerations are potentially significant for tracking 
the CO distribution near the photosphere because the limb is 
the region directly probed by transmissi on spectroscopy. We 
attempt to account for this by modifying Cooper & Showman 
(2005)'s equation ( 2). The new profile, shown in Figure|2] is 
identical to that of Cooper & Sho wmanl J2005I) over most of 
the day and night sides. But it has a smoother variation of 
temperature across the terminator, which remains fairly hot 
90° away from the substellar point (where half the stellar disk 
is still visible). It also increases the total radiative energy in 
the radiative-equilibrium system, which is realistic owing to 
the effects described above. 

5. SIMULATION RESULTS 

For this study, we have performed a suite of simulations 
with different input parameters, which all support the general 
conclusions of this work (^8}. For most of this discussion, we 
focus on the results of two simulations in particular. We call 
them the "nominal" and "cold" models, respectively. Other 
simulations are alluded to in this section and in ©|8| 

The nominal m odel contains input paramet ers identical 
to those used by ICooper &~ Showman (2005), which are 
tuned for HD 209458b specifically (although the radiative- 
equilibrium profil e ha s been slightly adjusted to account for 
limb effects; see ^4.4> . For the nominal model, we have used 
1000 K for the radiative-equilibrium temperature difference 
between the substellar point and the n ightside (at the top of 
the atmosphere, as described in SI4.3> . The cold model is a 
parameter variation on the nominal model. It is representa- 
tive of an imaginary EGP that is cooler than HD 209458b. 
The radiative-equilibrium temperatures for the cold model are 
similar to those expected for the transiting planet TrES-1, 
which orbits a K-type star (tFortnev et alJl2005t iBarman et alJ 



2005). 

We use for the cold model a new radiative-equilibrium pro- 
file, Tcoid(p)- This profile is 300 K cooler at every pressure 
level than the profile used for the nominal model: Tc \&(p) = 
Tim(p) - 300 K. For the cold case, we have also decreased 
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FIG. 3. — Temperature (grayscale) and winds (arrows) for our "nominal" 
model. This is a snapshot of the simulation at f s ; m = 1000 days. Maximum 
wind speeds in each frame are (a) 6600 m s , (b) 4300 m s , and (c) 3500 
m s . This is the HD 209458b nominal model, which uses parameters ap- 
propriate for the atmosphere of that planet. 

the assumed radiative-equilibrium temperature difference be- 
tween the substellar point and the nightside; for this model, 
we use Ar eq = 700 K. This is because we would expect cooler 
planets to have lower radiative forcing i n gen eral. With these 
adjustments, the procedure described in jj4.4l was rerun on the 
cold profile to determine the 3D radiative-equilibrium profile 
for the cold case. At the topmost layer (1 mbar), the substel- 
lar and nightside profiles are r ss = 1085 K and r n ; g ht = 385 K, 
respectively. 

We use the same values for the radiative relaxation time 
constant, T ra d(p), as for the nominal case. For a much 
cooler planet than HD 209458b, r, a d(p) would have to be re- 
evaluated, since radiative processes are temperature depen- 
dent. But this "cold" model is not enough colder than the 
nominal case for the radiative heating rates to differ by more 
than a factor of ~2. All other input parameters (acceleration 
of gravity, heat capacity, mean molecular weight, and radius) 
have been kept the same as in the nominal model. It should be 
emphasized, therefore, that the cold model is not a model of 
TrES-1 or any other known EGP in particular. It is simply a 
parameter variation on the nominal model, which is intended 
to elucidate the role of atmospheric temperature in determin- 
ing the CO/CH4 chemistry. 

For all simulation runs presented here, the horizontal grids 
contain 72 points in longitude and 45 points in latitude. We 
have run the models for f s j m = 1000 Earth days. This is long 
enough for the atmosphere model to equilibrate on layers hav- 
ing p < 3 bar, since the radiative time constants are fairly 
short — less than a month — in this region. We have verified 
that the kinetic energy in these low-pressure layers does not 
increase with time after 500-1000 days. We use a fairly small 
time step for all simulations: Af = 50 s. The small time step 
helps to control grid-scale instabilities. We discuss these is- 
sues further in 

We assume zero initial zonal and meridional winds: (u, v) 
= (0, 0). We also start the simulations with the temperature on 
every layer isothermal and equal to T\ m . We have verified that 
our solutions for the flow geometry — in the layers of interest 
here — hold even in the case of a n initially strong retrograde 
wind (Cooper & Showman 2005). For the initial tracer con- 
centration, we assume chemical equilibrium of CO at t slm = 0. 
We w ill discuss the effects of the initial tracer distribution in 

5.1. Temperature and Winds 

We focus here on results from the nominal and cold mod- 
els in the 10-1000 mbar range. These are the "visible layers" 
of EGP atmospheres; i.e., the layers that directly contribute 
to the outgoing radiation from the planet to space. The re- 
sults of our integrations for the nominal model can be seen 
in Figure|5] Figure [3] spans about 5 scale heights in pressure, 
encompassing the region crucial for observations of cl ose-in 
EGPs. Despite adjustments to the cooling scheme (see 34.4K 
the temperatures, wind speeds, and flow geometry of the nom- 
in al model are broadly consist ent with simulations published 
bv lCooper & Showman! (120051) . 

Strong temperature contrasts on isobars are evident in the 
simulation for all three layers shown. As seen in panel (a) of 



FIG. 4. — Nominal simulation, with inputs appropriate for HD 209458b. 
Shows the fraction of gaseous carbon present as CO in the (incorrect) equilib- 
rium picture of carbon chemistry. The other important carbon-bearing species 
is CH4 . Together, CO and CH4 comprise the total carbon budget of the at- 
mosphere iLodders & Feglev 2002). The top and middle panes show large 
gradients in the equilibrium concentration of CO, which correspond to steep 
horizontal temperature gradients (see Figure^. At 990 mbar, however, tem- 
peratures are everywhere higher than 1200 K. On most of this layer (bottom 
pane), CO is therefore thermodynamically favored over CH4 . 



FIG. 5. — The fraction of carbon as CO in the nominal model, in which we 
have coupled the chemical kinetics to the dynamics. The simulation results 
emphasize that chemical equilibrium does not hold in this region of the atmo- 
sphere. The simulation predicts a homogeneous distribution of CO, with CO 
concentrations in local cool regions exceeding chemical equilibrium values 
by many orders of magnitude. 

Figure |3 the atmosphere is nearly in radiative equilibrium at 
low pressures, where the radiative time constants are less than 
a day. At higher pressures, panels (b) and (c), the radiative 
heating is weaker and large departures from radiative equilib- 
rium can be seen. A strong superrotating jet extending from 
the equator to the mid-latitudes develops in the flow. In the 
middle panel, (b), the hottest regions of the atmosphere are 
blown downstream from the substellar point (A = 0°, <j> = 0°) 
by about 45°. This layer at 106 mbar is near the planet's 
photosphere. Longitudinal temperature contrasts are appar- 
ent near 1 bar, as shown in panel (c), where the hot region of 
the atmosphere is shifted downwind from the substellar point 
by a full 90°. 

Many features of the flow in the cold model (not shown) are 
quite similar to the nominal model, though the temperatures 
in this atmosphere are on corresponding isobars are notably 
cooler (by ~ 300 K). From a dynamics standpoint, the models 
are basically the same, with similar flow patterns and wind 
velocities developing in each case. The temperature differ- 
ences between the nominal and cold models are important for 
determining the carbon chemistry, however, since the reaction 
kinetics depend sensitively on temperature (see ^3}. 

5.2. CO Distribution 

Our result for the distribution of gaseous carbon in HD 
209458b's atmosphere is shown in Figure |5] Contrast this 
with Figure^] which shows how carbon would be distributed 
under (fictitious) chemical equilibrium conditions. Grayscale 
in these figures represents the percentage of total atmospheric 
carbon present as CO. 

As shown in panels (a) and (b) of Figure 0] the assump- 
tion of chemical equilibrium implies that vast regions of HD 
209458b's atmosphere are devoid of CO. In these cool re- 
gions, CH4 is strongly favored thermodynamically. Steep gra- 
dients in the concentration of CO are possible wherever the 
temperature varies rapidly ab out the line of equal CO /CH4 
abundance (refer to Figure 2 of Lodders & Feglev (2002)). By 
contrast, panel (c) shows that near the 1 bar level in the atmo- 
sphere, CO is the more thermodynamically stable compound. 
This is because near 1 bar, as shown in Figure [3] the tem- 
perature is everywhere higher than 1200 K, which is where 
the equilibrium abundances of CO and CH4 are equal (at this 
pressure). 

Our model's result, which considers the extremely slow rate 
of CO / CH4 interconversion in this region of the planet's at- 
mosphere, shows a completely different view of the carbon 
chemistry. The simulation predicts a nearly homogeneous dis- 
tribution of CO in all three layers shown. CO permeates the 
atmosphere in this region, being by far the dominant carbon- 



Carbon Chemistry on HD 209458b 



9 



bearing species (at the 98-99% level). 

Our result for the cold model is similar to the nominal case. 
The cold atmosphere simulation also predicts that carbon is 
homogenized in the upper atmosphere. As in the nominal 
case, carbon is the major carbon-bearing gas in the 10-1000 
mbar range. Here, though, 20% of the carbon is CH4 (as op- 
posed to 1-2% in the nominal simulation). The figure reveals 
that CO is likely to be highly abundant near the photospheres 
of EGPs — both on the dayside and nightside — even on planets 
like TrES-1 that are significantly cooler than HD 209458b. 

We have also run a "hot model," in which 300 K was added 
to the radiative-equilibrium profile of the nominal model: 
THot = Ti ro + 300 K. The circulation pattern seen for the nom- 
inal and cold models (Figure[3} develops for the hot model as 
well, though the temperatures on each corresponding layer are 
higher by ~ 200-300 K. The temperature is high enough in the 
hot model for CO to be present in high abundance everywhere 
in both the equilibrium and disequilibrium pictures. There- 
fore, planets significantly warmer than HD 209458b should 
contain no CH4 at all near their photospheres. 

Steady-state has been reached for p < 3 bar. But the ki- 
netic energy of the deeper layers in the model still increases 
with time after 1000 days. Large equator-to-pole temperature 
variations are seen in the bottom three panels (spanning the 
3-10 bar range). Wind speeds remain high at these pressures. 
Although heating has been turned off for pressures exceeding 
10 bar, the atmosphere is not motionless even at 30 bar. We 
attribute this to vertical transport processes that transfer ki- 
netic energy and angu lar momentum to the p > 10 bar region 
(Cooper & Showman 2005). 

The vertical wind in the model is quite strong (u> = 0.05- 
0. 1 mb s" 1 at 3 bar). The vertical wind shows evidence of both 
small-scale (~ 10°) and large-scale (~ 45°) structure, with al- 
ternating regions of upwelling and downwelling. The small- 
scale features are likely the result of atmospheric waves. The 
large-scale structures correspond to overturning circulations 
superimposed on the fast superrotation, although we empha- 
size the region is not convective (vertical motion on this planet 
is physically much more similar to stratospheric overturn- 
ing in the Earth's atmosphere). Upwelling occurs primarily 
within 20° of the equator, whereas prominent downwelling 
occurs in the mid-latitudes. In the 1-1000 mbar region, there 
is also a longitudinal dependence of the vertical wind. In these 
upper layers of the model, strong upwelling occurs in hot re- 
gions (which is downstream of the substellar point at 1 bar), 
whereas downwelling occurs in cool regions. This longitu- 
dinal structure is present but much weaker in the 1-10 bar 
region. 

6. ANALYSIS / INTERPRETATION 

Ultimately, the high CO concentrations we see in these 
simulations are the result of chemical quenching. In gen- 
eral, the distribution of carbon-bearing species will depend on 
(1) where in the atmosphere quenching occurs (the "quench 
level"), and (2) the chemical equilibrium abundances of CO 
and CH4 at the quench level. The strong vertical wind is ca- 
pable of transporting parcels through a pressure scale height 
(~ 500 km) in ~ 10 5 s, which is comparable to the timescale 
for horizontal advection of parcels across a hemisphere. This 
supports the notion that quenching occurs primarily in the 
vertical, not the horizontal, due to uniformly long chemical 
timescales in the region of disequilibrium (see §0. 

The quench level of this model — which demarcates the re- 
gion where the equilibrium and disequilibrium views give the 



same answer — is in the p = 3-5 bar range. This is of course 
not exact because quenching depends on temperature, and 
considerable equator-to -pole temperature va riations exist at 
depth in these models (Cooper & Showman 2005). For ex- 
ample, comparison of the bottom panels of Figures |4] and [5] 
reveals that at the equator, quenching likely occurs closer to 
p = 2 bar. But at higher latitudes, where the temperature is 
lower on these layers, the quench level is deeper {p = 6 bar). 

The cause of such high CO abundances from 10-100 mbar 
becomes clear. The value of Xco in equilibrium is quite high 
in the region of quenching. It is in fact generally higher near 
3 bar than it is in deeper layers (not shown). This is due to the 
thermodynamic stability of CH4 at high pressures, as shown in 
Table [2 The results clearly indicate, though, that the quench 
level (not deeper layers) determines the photospheric value. 
This is due to vertical transport. The strongest vertical wind 
is near the equator, where the temperature is also the highest 
(on a given vertical layer). It is for this reason that the dis- 
equilibrium values of CO in the 1-1000 mbar region are so 
uniformly high: vertical mixing is maximum at the equator at 
the quench level, where temperature and pressure conditions 
strongly favor the formation of CO. 

We now investigate a simpler view of quenching, which 
yields very similar results to the numerical simulations. 

6.1. Eddy Diffusion 

The simplest model of disequilibrium chemistry defines the 
quench level to be where the timescale for reactions equals the 
timescale for dynamical advection in the system: i.e., quench- 
ing occurs at T crie m = Tdyn- The latter quantity, Td yn , represents 
the timescale for dynamical motions to advect material across 
regions of atmosphere differing appreciably in temperature 
and pressure. Depending on local gradients in temperature, 
these timescales can be either large or small. 

These considerations led Smith JT9981) to revise the 
eddy diffusion model o f quenching originally proposed by 
iPrinn & Barshay] (119771) . In ISmithi (Il998j) 's model, the dy- 
namical timescale is determined by eddy diffusion: 

Tdyn = — , (23) 
W 

where L is a length scale and w is the vertical velocity near 
the quench level. 

Previous researchers estimated the quench level on Jupiter 
assuming L = H, where H is the pressure scale height of the 
atmosphere. Smith ( 1998) demonstrates that the assumption 
of L = H is incorrect for most planetary atmospheres. How- 
ever, the eddy diffusion calculation is approximately correct 
if L is replaced with an effective length scale L e g, "which de- 
pends on the e-folding length scales of r C hem and Td yn and the 
pressure dependence of the equilibrium value of the chemical 
abundance." He goes on to develop a recipe for determining 

L e g- 

We apply the ISmTthl ( fT998i) recipe, which offers excellent 
insight into the chemistry governing our simulations, to the 
nominal atmosphere model. We also assume for the moment 
that quenching occurs in the vertical direction, not horizon- 
tally. The crucial quantity needed to estimate eddy diffu- 
sion is the vertical velocity (in height coordinates). Given 
u> = Dp/Dt, the vertical velocity in height coordinates is ap- 
proximately given by 

w = u/(-pg), (24) 

where p is the density and g is the acceleration of gravity. 
It should be noted that this is not exact but is an excellent 
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FIG. 6. — De monstrates quenching according to the eddy diffusion model 
of Smith i 1998) in the nominal simulation for two vertical columns of atmo- 
sphere. The top panel shows T c hem and T& yri on two vertical columns: (1) 
(A = 120°, = 0°) (solid lines) and (2) (A = -60°, = 0°) (dashed lines). 
The quench level for the first column is where the thick and thin solid lines 
cross. Likewise, the quench level at (A = —60° , <fi = 0°) is where the thick and 
thin dashed lines cross. The middle and bottom panes show the temperatures 
and vertical velocities used to calculate rjy,, in each region. 



FIG. 7. — Comparison of actual (solid line) and equilibrium (dashed line) 
CO abundance in the nominal simulation. Plotted is the % of carbon as CO — 
as opposed to CH4. Titles above the panels give the longitude and latitude for 
each column. The horizontal dotted lines show the quench level for each 
column as predicted with the Smith 1 1998) recipe. 



approximation for shallow atmospheres jHoltonl 19921 p. 77- 
80). 

To estimate quenching, it is essential to determine the 
appropriate length scale L e g to use in comp uting Tdm, the 
timescale for eddy diffusion. Following the ISmithl dl998l) 
recipe, L e ff ~ 0.5 H in the first region, and Leg ~ 0.3 H second 
region. These are closer to H than are appropriate for Jupiter 
because the equilibrium abundance of CO is uniformly high 
at the equator for p > 200 mbar. This correction nonetheless 
causes a significant shift in the pressure at which quenching 
occurs in each case. 

We show the quench levels (using Td yn = L e ff/w, not H/w) 
for both profiles in the top panel of Figure[6] The quench level 
near (A = 120°, (j) = 0°) is about 3.0 bars; near (A = -60°, = 
0°), it is at 2.5 bar. The middle and bottom panels of Figure 
[6] show the vertical velocities and temperatures in these two 
regions. Here, we have averaged columns of w over a 45° 
angular radius on the sphere, which removes small-scale vari- 
ations but preserves the large-scale structure of the vertical 
wind. For p < 1 bar, downwelling occurs in the first region, 
whereas the second region, which is centered on the opposite 
hemis phere, has upw eiling. 

The Smith ( 1998) quenching model also predicts a homo- 
geneous distribution of CO above the quench level, with the 
value determined by the equilibrium CO abundance at the 
quench level. Since this is high in both cases, the eddy dif- 
fusion model also predicts very little CH4 at the photosphere 
of HP 209 458b. To emphasize the agreement between the 
ISmitr]lfT998l) view of quenching and the actual abundance of 
CO in our simulations, we compare the equilibrium view to 
the CO abundance predicted in our simulation in Figure [7] 
This has been done on four columns of atmosphere to show 
the variation in quenching in different atmospheric regions. 
The Smith ( 1998) quench level in each pane is indicated by 
a thin dotted line, which agrees well with the convergence of 
the equilibrium vs. disequilibrium profiles (thick dashed and 
solid lines, respectively). 
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We conclude this discussion by reiterating why horizontal 
quenching near the day-night terminator is not of great im- 
portance in the simulation. Using reasoning analogous to our 
discussion above, the horizontal mixing timescale should be 
Tdyn.h = Lh/u, where // is the supersonic zonal wind and Lh is 
a length scale appropriate for horizontal quenching. Lacking 
better information, Lh = 2irR p (just the circumference of the 
planet). Note that the true value of the horizontal mixing time 
may be smaller than 2ttR p /u in regions where parcels travel 
through large temperature gradients. But given the extremely 
fast wind speeds (u ~ 3000 m s" 1 ) predicted in the upper lay- 
ers, Tdyn.h can never be larger than this value. Assuming for 
the moment that this is correct, we obtain Td yn ,h ~2x 10 5 s. 
But Figure |6]clearly shows that r c h e m 3> 10 s for p < 3 bar, 
where vertical quenching is achieved. Hence, the quench 
level in the vertical has already been reached at levels where 

7"chem 

7. DISCUSSION 

7.1. Role of Photochemistry 

iLiang et all (12004 ) show that the distribution of photochem- 
ical products, including aerosols and radicals, is insignificant 
at the pressure levels considered here compared to the primary 
constituents (H2, He, H2O, CO, and possibly CH4). But the 
effects of these species on the CO/CH4 reactions has not been 
explored here. Conceivably, OH and H coming from the pho- 
tolysis of H2O (Liana et al. 2003) could influence the reaction 
rates of the Yung et al. ( 1988) sequence. 

We have not explored the consequences of photochemically 
enhanced OH and H populations on the CO / CH4 interconver- 
sion chemistry here. This effect is most likely to be notice- 
able in the highest layers of the model (p < 10 mbar), above 
which the stellar ultra-violet (UV) light is absorbed. In or- 
der to have an important effect on our predictions for the CO 
distribution, however, the chemical interconversion timescale 
resulting from interaction with photochemical radicals would 
have to be reduced to r C h e m = 10 6 s (or less). Such timescales 
would be very short compared with the values listed in Table 
|2]at these temperatures and pressures. That is, unless the pho- 
tochemical production of radicals dominates the chemistry, 
the disequilibrium distributions will be substantially similar 
to those presented here. 

7.2. Treatment of Radiation 

Asynchronous rotation on hot Jupiters is possible, either 
due to eccentricity/obliquity pumping by other plane ts in 
the system or thermal tides ( Showman & Guillotl 120021) . In 
the case of HD 209458b, we expect synchronous rotation to 
hold more closely than planets in eccentric orbits. An asyn- 
chronous component to the rotation, due to thermal tides and 
transfer of angular momentum from the atmosphere to the in- 
terior, could conceivably influence the atmospheric dynamics 
solution (Figure [3}- For example, our time-steady flow sce- 
nario might be superposed in some regions by a slow period- 
icity resulting from a radiative heating profile that slowly ad- 
vances along the meridians. But as we discuss in |6] quench- 
ing is in the vertical. Adding asynchronous rotation to the 
model might change the details of the wind and tempera- 
ture patterns but not the mean vertical velocity. Hence, the 
CO/CH4 distributions should not be sensitive to this correc- 
tion. 

Newtonian heating ignores potentially subtle feedbacks in 
the radiation physics. It has been used successfully, how- 



ever, to model essential features of Earth's general circulation 
llHoltonll 19921 p. 409-411). We note that the approximation 
is likely to be less valid deep in the atmosphere of EGPs than 
it is for Earth. This is due to long radiative time constants 
and huge departures from radiative equilibrium. Hence, the 
use of Newtonian heating is probably the major source of er- 
ror in our atmosphere model. We use it here partly because 
the approximation is reasonably valid for p < 2 bar, where 
T ra d is less than several days. The scheme is also relatively 
easy to implement and performs well computationally, which 
is very desirable, as the AGDC2 itself demands considerable 
computational resources. We are able, therefore, to test the 
sensitivity of the results to a broad set of input parameters. 

We n ote also that th e abundances of H2O, CO, and CH4 — 
used by Ir o et al. I J20051 to derive the contributions of these 
molecules to the total opacity — are calculated based on chem- 
ical equilibrium models. We demonstrate here, however, that 
chemical equilibrium does not hold for p <C 10 bars. This in- 
ternal inconsistency is inherent in our use of Newtonian cool- 
ing and is unfortunately unavoidable. Radiative-equilibrium 
models of HD 209458b's atmosphere to date have all assumed 
chemical equilibrium. Even in the future, this coupling be- 
tween dynamics and radiation will be difficult to disentan- 
gle because, as has been previously discussed, disequilibrium 
a bundances depen d on the overall meteorology. 

Ilro et al.l d2005l> also do not consider clouds. Where present, 
high silicate clouds could dominate the opacity and greatly 
perturb the atmospheric profile. Cloud deck optical thick- 
ness in general varies with wavelength, depending on the 
vertical extent of the cloud and the condensate particle sizes 
( Ackerman & Marlev 2001; Coope ret alJ2003l) . It is unclear 
what properties a silicate cloud high in the atmosphere would 
have. This is a complicated subject involving intricate feed- 
backs between radiation, cloud microphysics, and dynamics. 
We neglect the possibly crucial effects of clouds here. Includ- 
ing condensation processes will be a major future advance in 
the field of EGP atmosphere modeling. 

7.3. Initial Conditions 

Given the large r c h em values at low pressures, it is rea- 
sonable to wonder whether our solutions for the CO dis- 
tribution depend on the initial conditions. For the simula- 
tions shown, we have assumed that CO is initially in chem- 
ical equilibrium everywher e (equation As discussed in 
Cooper & Showman (2005), temperatures and winds are not 
strongly sensitive to the initial conditions. The primary con- 
cern here is whether vertical transport can mix CO upward 
rapidly enough to bring the CO mole fraction at the photo- 
sphere up to its value at the quench level. 

Clearly, chemical processes occurring on timescales longer 
than ~ 1 year cannot have run to completion by 1000 days 
of integration time. But we can argue on simple grounds 
that 1000 days is still long enough to determine Xqo, even 
where r c h em 3> 1 year. Assuming the vertical wind at 1 bar 
is ~10ms _1 , as shown by the model, and the vertical depth 
from the photosphere (100 mbar) to the quench level (3 bar) 
is ~ 2% of the planetary radius (several scale heights), we 
would expect upward mixing to occur on timescales of ~ 1 
week. Thus, after ~ 1000 days of integration time, the CO 
concentration should approach the value at the quench level. 

We have run two additional simulations with different ini- 
tial conditions. The first simulation assumes no initial CO in 
the atmosphere; i.e., 100% of the carbon is CH4. The second 
simulation assumes no initial CH4; i.e., 100% of the carbon 
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is CO. Both use the input parameters of the nominal model in 
all other respects. Therefore, temperatures and winds are the 
same as in the nominal case. After even several hundred days 
of simulation time, the solutions for CO in these two cases are 
virtually indistinguishable from the nominal simulation (see 
Figure|5}. These tests demonstrate that memory of the initial 
condition for Xqo has been erased by f s ; m = 1000 days. 

7.4. Metallicity and C/O 

Our nominal s imulation assume s solar abundances based 
on the Lodders & Fealev ( 1998) compilation. We note 
that the precise values of the solar abundances, especially 
C, N, and O, have undergone minor modifications in the 
past seven years, some of which are still under discus- 
sion within the scientific community dGrevesse & Sauval 

e 98t iMahaffv et al.ll2000t lAllende Prieto et alJl2001t 120021 
iddersl l2003l lAsplund et al.ll2006l) . In particular, the CNO 
abundances, along with those of the noble gases, have been 
revised down from the ge o chemically based a bundances of 
lAnders & Grevessd (U989T) : fodders & Fegley ( 1998). The 
C/O ratio is now thought to be exactly 0.5; following 
lLodders & Feelevl 0998), we use C/O = 0.57 here. Incor- 
porating these updates to the solar abundances will modify 
the results presented here somewhat, though the correction is 
likely to be fairly minor. 

From stellar spectroscopy, the endowment of heavy ele- 
ments in the star HD 209458 i s well-known to be about solar: 
[Fe/H] = 0.014 ±0.03 jFischer & Valentill2005l) . That says 
nothing about the metallicity of the planet orbiting it, which 
has not yet been measured. Jupiter and Saturn boast pro- 
nounced enrichments in heavy elements relative to the Sun, 
suggesting EGPs may also be metal-rich compared to their 
host stars. We consider here the effects of increasing the rel- 
ative proportions of carbon and oxygen in HD 209458b's at- 
mosphere. 

Suppose the constants c\ and C2 in equations ^ an d ID are 
both increased by a factor of 10. This is a much greater metal 
abundance than is seen on Jupiter, which has enrichment of 
carbon and nitrogen by a factor of ~ 3 relative to the Sun (from 
Galileo probe measurements; the probe w as not successful in 
measuring the abundance of oxygen — see IWons et al]l2004l) . 
As carbon and oxygen in this case are still minor constituents, 
this change will not appreciably affect the mean molecular 
weight or heat capacity of the atmosphere. According to our 
chemical equilibrium solution, equation |5] increased [Fe/H] 
favors the formation of CO. That is, for a given temperature 
and pressure, the ratio Xqo/Xch 4 is highe r for r Fe/Hl = +1.0 
than it is at [Fe/H] = 0.0, in agreement with Lodders & Feglev 
( 2002). This implies that, relative to the solar metallicity case, 
the quench level should have a greater abundance of CO in 
the [Fe/H] = +1.0 atmosphere by more than a factor of 10. 
Hence, we would also expect the photosphere to show highly 
abundant CO. 

We have run a high-metallicity simulation, in which we set 
[Fe/H] = 1.0 in the atmosphere. The other parameters of this 
model are identical to the nominal simulation, so the dynam- 
ics are the same. We indeed see in this simulation that 100% 
of the carbon atoms are in CO, with no CH4 at all (as opposed 
to 1-2% CH4 predicted in the nominal case). 

We stress that the above discussion assumes atmospheric 
temperatures do not change drastically as a result of increas- 
ing the metallicity. This assumptio n may not be corr ect for all 
EGPs. We note, however, that Fortnev et al. (2005) consider 
metallicity adjustments to their radiative-equilibrium models 



of HD 209458b and TrES-1. They find that the effective tem- 
perature of TrES-1 increases by only ~ 14 K for a metallic- 
ity of 5x solar. The direction of the effect — warming, not 
cooling — also favors the appearance of CO. 

Similarly, increasing the C/O ratio in the atmosphere 
will not inhibit CO's formation at the temperature of HD 
209458b's atmosphere. In the Sun, C/O w 0.5. If C/O > 1, 
CO can potentially take up 100% of the available oxygen, 
in which case H2O will be depleted. Atmospheres showing 
the spectral signatures of CO without prominent H2O features 
very likely have a high C/O ratio relative to solar. The effects 
of super-solar abundances and enhanced C/O are discussed 
in mo re detail elsewhere ( Fortne vet alJ 120051: ISeager et alJ 
120051) . 

8. CONCLUSIONS / IMPLICATIONS 

We have run 3D numerical simulations of HD 2094548b's 
atmosphere to explore the chemistry of carbon. We advect 
CO through the atmosphere as a passive tracer, relaxing its 
concentration toward the local chemical equilibrium value. 
The relaxation timescale depends on the reaction pathways 
that convert CO to CH4. We follo w recent previous work 
on disequilibrium carbon chemistry ( Griffith & Yelle] | 19991: 
Bezard et all2002l) and assume that the Yung et al. ( 1988) se- 
quence is the most efficient mechanism for the chemical re- 
duction of CO. 

We calculate the relative abundances of CO, CH4, and H2O, 
which are the main species of carbon and oxygen in the atmo- 
sphere. Our nominal simulation models HD 209458b specifi- 
cally. For comparison, we have also run simulations of "cold" 
and "hot" atmospheres. The cold simulation is representative 
of a planet that is ~ 300 K cooler than HD 209458b; the hot 
simulation models a planet ~ 300 K warmer. 

The atmospheric dynamics by 1000 Earth days of sim- 
ulatio n time has reached the steady state described in 
Cooper & Showman (2005). Eastward winds in the nomi- 
nal model reach ^4000 m s" 1 at the photosphere, with corre- 
sponding temperature contrasts of ~ 500 K. Atmospheric dy- 
namics partially redistribute the thermal energy to the planet's 
nightside, but longitudinal contrasts in temperature are appar- 
ent down to ^2 bars. 

Chemical equilibrium has been the basis of all prior dis- 
cussions about CO and CH4 on HD 209458b. We show here, 
though, that the equilibrium picture of CO/CH4 chemistry is 
a vast misconception. In regions where chemical equilibrium 
does not hold, the true CO abundance depends not only on 
local quantities but also on the overall meteorology. 

We show that the timescale for CO/CH4 interconversion 
depends sensitively on both temperature and pressure (Table 
|2ji. At low temperatures and pressures, reactions are slow; 
at high temperatures and pressures, reactions are fast. As a 
result of long reaction times, CO is not expected to be in 
chemical equilibrium in the 1-1000 mbar region (the "visi- 
ble" layers). Chemical equilibrium in the model atmosphere 
is attained deeper in the atmosphere (p > 5 bars), where the 
interconversion timescale is short relative to dynamical mix- 
ing times. 

The atmosphere is well-mixed by vertical motions, which 
leads to vertical (as opposed to horizontal) quenching. Verti- 
cal quenching results in photospheric CO concentrations ex- 
ceeding chemical equilibrium values by many orders of mag- 
nitude. The simulations show that CO, CH4, and H2O will 
be homogeneously distributed throughout the photosphere of 
HD 209458b and similar planets. At p < 2 bars, almost all 
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of the carbon is in the form of CO, even in cold regions such 
as the nightside. Consequently, we expect water vapor to be 
less plentiful than previously thought. These conclusions are 
consis t ent w ith the simpler view of quenching presented by 
Smith! (Tf998l) . 

The metallicity of HD 209458b is not known. Our nomi- 
nal simulation assumes solar abundances. We have also run 
a variation on the nominal model that assumes ten times en- 
hancement of the heavy elements: [Fe/H] = +1.0. Disequilib- 
rium effects operate similarly in this model, and the ratio of 
CO to total carbon in the atmosphere remains close to 100%. 
Indeed, higher metallicity favors the formation of CO over 
CH4. Hence, if close-in EGPs have characteristically super- 
solar abundances, CO could be prominent at the photospheres 
of planets even cooler than T e ff = 1000 K. 

Likewise, increasing the C/O ratio would not significantly 
inhibit CO's formation at depth (and subsequent transport up 
to the photosphere). If C/O exceeds the solar value of 0.5 
in the Sun, CO on HD 209458b should be highly abundant. 
The tell-tale sign of this scenario at HD 209458b's effective 
temperature is H2O depletion. 

Our results have several interesting consequences for near- 
future observations of close-in EGPs, although a thorough 
analysis of close-in EGP spectra is beyond the scope of this 
work. We conclude simply with a brief description of the im- 
plications of these findings. 

The hot model shows that planets significantly warmer than 
HD 209458b should have no CH4 in their photospheres, un- 
less the C/O ratio is greater than 1. But in planets signifi- 
cantly cooler than HD 209458b, disequilibrium effects should 
raise the concentration of CO to detectable levels. The spe- 
cific value of Xqo cannot be predicted with definite precision 
based on these simulations; there are too many uncertainties 
in the model's inputs. However, examination of the values of 
fchem (Table |3 suggests that 1-10 bar may be the approx- 
imate quench level for all close-in EGPs within 750 K of 
HD 209458b's effective temperature (T eff ~ 1350 K for this 
planet). The relative concentrations of CO and CH4 at the 
photosphere are therefore diagnostic of local conditions in the 
1-10 bar region of the atmosphere, which is deeper than ob- 
servations can directly probe. Neith er species has y et been de- 
tected in this planet's atmosphere (Richards on et alJ l2003b a; 
iDeming et alJl2005al) . Our results here suggest that CH4 de- 
pletion (not depletion of CO) is the more likely situation; even 
in chemical equilibrium, CH4 is likely to be depleted at the 
limb. 

Measurements of carbon-bearing species can thus serve as 
a useful calibration tool of radiation and dynamics models of 
these planets. It should be noted that metallicity also plays a 
role in determining the abundance of CO at the quench level. 
This may be a difficult effect to deconvolve from tempera- 
ture and pressure, unless absolute concentrations can be mea- 
sured. For near-future explorations, what is crucial here is that 
planets with little or no CH4 must be very much warmer than 
TrES-1 or else cooler, but very metal-rich (as high metallicity 
favors the formation of CO over CH4). 

It is known that HD 209458b do es contain carbon and oxy- 
gen JVidal-Madiar et alJ20"04h . As lDeming et alJ ll2005al) dis- 
cuss, the null detection of CO using NIRSPEC on Keck II 
cannot be easily attributed to CO depletion in HD 209458b's 
atmosphere. This lends credence to the hypothesis of a high 
cloud or haze obscur ing the spectral features of CO near 2 
/im. As discussed by Fortnev (2005), however, the slant ge- 
ometry implies that even a layer of modest optical thickness 



could hide the signature of CO. It may be possible to distin- 
guish between these and other scenarios. 

Our results suggest that TrES-1, which orbits a K-type star, 
should also show strong CO absorption features in its spec- 
trum. A high silicate condensation cloud obscuring the pho- 
tosphere is not expected on TrES-1 due to lower temperatures 
overall. The condensation of silicates should occur at much 
greater depth in TrES-1 's atmosphere. It would be very illu- 
minating to perform a transmission spectroscopy search for 
CO on TrES-1. Discovery of CO on TrES-1 would make a 
strong case for high silicate clouds in the atmosphere of HD 
209458b. If CO is also not detectable on TrES-1 by transmis- 
sion spectroscopy, on the other hand, the hypothesis of strato- 
spheric haze (that hides CO on both planets) becomes more 
appealing. 

The ubiquitousness of CO throughout the photospheres of 
close-in EGPs should also have profound implications for the 
spectra of these objects. Our results for the cold model sug- 
gest that CO may even be abundant on cooler EGPs than 
TrES-1. If so, CO's spectral signature should be detectable. 
The absorption features of H2O may also be less pronounced, 
especially on hot planets, in which over half of the avail- 
able oxygen atoms are boun d up in CO. In addi tion to the 
2 /im band pass explored by Deming et al. (2005a), we sug- 
gest searches for CO in the rovibrational fundamental near 4.5 
micron. If present, the opacity of CO overwhelms that of H2O 
vapor at this wave length, even more so than the wavelength 
region explored in Demi n~et alJ l l2005al) . As an IR window 
also exists at 4.5 /im, this observation may also be possible 
using ground-based platforms. 

It is noteworthy also that these simulations predict a homo- 
geneous distribution of CO, CH 4 , and H 2 in HD 209458b's 
atmosphere, although the tempe rature is almost certainly not 
isothermal at the photosphere. Coop er*& Showmai] (2005) 
show that temperature variations over the photosphere can 
lead to brightness variations in the planet with orbital phase. 
It must be emphasized t hat the synthetic light curve of 
Cooper & Showman ( 120051) pertained to the bolometric lu- 
minosity of the planet (i.e., the emission integrated over all 
wavelengths); nothing is directly implied about the planet's 
brightness at a particular wavelength. The flux contrast of the 
planet's photosphere over a full orbital period would likely 
be suppressed somewhat by a homogeneous distribution of 
CO and CH4, as the distribution of these species determines 
the level of the photosphere [J. Fortney, private communica- 
tion]. Wavelength-dependent light curves will require gener- 
ating synthetic spectra using a radiative transfer model, a task 
we leave for future work. But we hope this discussion will 
inspire near-future searches for CO and CH4 on these planets, 
some of which should be possible even with ground-based 
telescopes. 
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Carbon Chemistry on HD 209458b 



TABLE 1 
CO in Chemical Equilibrium 



T [K] 10 mbar 100 mbar 500 mbar 1 bar 5 bar 10 bar 100 bar 500 bar 1 kbar 



500 


1.7 x 10~ 9 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 




0.0 




0.0 




750 


5.5 x 10~ 2 


6.0 x 10~ 


4 2.4 x 10~ 5 


6.0 X 10~ s 


2.4 x 10~ 7 


6.o x ur 8 


6.3 x 10" 


-10 


0.0 




0.0 




1000 


1.0 


0.73 


0.15 


4.4 X 10~ 2 


1.9 x 10- 3 


4.6 x 10- 4 


4.7 x 10- 


-6 


1.9 x 10- 


-7 


4.7 X 


io- 8 


1250 


1.0 


1.0 


0.95 


0.85 


0.27 


9.4 x 10~ 2 


1.1 x 10- 


-3 


4.4 x 10- 


-5 


1.1 X 


10- 5 


1500 


1.0 


1.0 


1.0 


0.99 


0.89 


0.71 


4.0 x 10- 


-2 


1.7 x 10- 


-3 


4.2 X 


10" 4 


1750 


1.0 


1.0 


1.0 


1.0 


0.99 


0.96 


0.32 




2.2 x 10- 


-2 


5.7 X 


IO" 3 


2000 


1.0 


1.0 


1.0 


1.0 


1.0 


0.99 


0.70 




0.13 




3.7 X 


ur 2 


2500 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


0.96 




0.61 




0.32 





NOTE. — Fraction of carbon in the atmosphere present as CO (i.e., the ratio Xco / c\ ), assuming chemical equilibrium (see eq.lTI— [5|. We obtain these 
values by solving equation|5j which calculates the chemical equilibrium of CO in a gas of solar composition. Our treatment follows Lodders & Feslev 
2002). Data for the standard Gibbs free energies of formation for CO, CH4, and H2O are taken from the NIST-JANAF tabulations Ichasen998ft . 



TABLE 2 

CO/CH4 INTERCONVERSION TIMESCALE 



T[K] 10 mbar 100 mbar 500 mbar 1 bar 5 bar 10 bar 100 bar 500 bar 1 kbar 



500 


1.0 x 10 24 


4.3 x 10 21 


1.7 x 10 20 


4.9 x 10 19 


3.7 X 10 18 


1.3 x 


10 18 


4.0 X 10 16 


3.5 x 10 15 


1.3 x 10 15 


750 


9.0 x 10 17 


2.8 x 10 15 


5.2 x 10 13 


9.3 x 10 12 


1.9 x 10" 


3.9 x 


10 10 


4.2 X 10 8 


3.1 x 10 7 


1.1 x 10 7 


1000 


1.1 x 10 15 


3.6 x 10 12 


6.4 x 10 10 


1.1 x 10 10 


2.1 x 10 8 


3.7 x 


LO 7 


1.4 X 10 s 


4.8 x 10 3 


1.3 x 10 3 


1250 


2.3 x 10 13 


7.4 x 10 10 


1.3 x 10 9 


2.3 x 10 s 


4.2 X 10 6 


7.4 x 


10 5 


2.4 x 10 3 


52 


11 


1500 


1.9 x 10 12 


6.0 x 10 9 


1.1 x 10 s 


1.9 x 10 7 


3.4 x 10 5 


6.0 x 


1() 4 


1.9 x 10 2 


3.6 


0.68 


1750 


3.3 x 10" 


1.0 x 10 9 


1.9 x 10 7 


3.3 x 10 6 


5.9 x 10 4 


1.0 x 


LO 4 


33 


0.60 


0.11 


2000 


9.1 x 10 10 


2.9 x 10 s 


5.1 x 10 6 


9.1 x 10 5 


1.6 x 10 4 


2.9 x 


10 3 


9.1 


0.16 


2.9 x 10~ 2 


2500 


1.6 x 10 10 


5.0 x 10 7 


8.9 x 10 s 


1.6 x 10 5 


2.8 x 10 3 


5.0 x 


10 2 


1.6 


2.8 X 10~ 2 


5.0 x 10~ 3 



NOTE. — Timescale in seconds for CO/CH4 interconversion in the! Yun g et a L, ■ 198 8) catalyzed reaction sequence at the same temper ature s and 
pressures used in Table [71 These values have been computed using equations n9l and |l8l following equations (7) and (8) of Bezard et al. 200.3) ■ The 
expression for r C hcm derives from consideration of the rate of the reverse of reaction|9| which has been measured in the laboratory 
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